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ABSTRACT 


Context.  Helicity  is  a  fundamental  property  of  magnetic  fields,  conserved  in  ideal  MHD.  In  flux  rope  topology,  it  consists  of  twist 
and  writhe  helicity.  Despite  the  common  occurrence  of  helical  structures  in  the  solar  atmosphere,  little  is  known  about  how  their  shape 
relates  to  the  writhe,  which  fraction  of  helicity  is  contained  in  writhe,  and  how  much  helicity  is  exchanged  between  twist  and  writhe 
when  they  erupt. 

Aims.  Here  we  perform  a  quantitative  investigation  of  these  questions  relevant  for  coronal  flux  ropes. 

Methods.  The  decomposition  of  the  writhe  of  a  curve  into  local  and  nonlocal  components  greatly  facilitates  its  computation.  We  use  it 
to  study  the  relation  between  writhe  and  projected  S  shape  of  helical  curves  and  to  measure  writhe  and  twist  in  numerical  simulations 
of  flux  rope  instabilities.  The  results  are  discussed  with  regard  to  filament  eruptions  and  coronal  mass  ejections  (CMEs). 

Results.  (1)  We  demonstrate  that  the  relation  between  writhe  and  projected  S  shape  is  not  unique  in  principle,  but  that  the  ambiguity 
does  not  affect  low-lying  structures,  thus  supporting  the  established  empirical  rule  which  associates  stable  forward  (reverse)  S  shaped 
structures  low  in  the  corona  with  positive  (negative)  helicity.  (2)  Kink-unstable  erupting  flux  ropes  are  found  to  transform  a  far  smaller 
fraction  of  their  twist  helicity  into  writhe  helicity  than  often  assumed.  (3)  Confined  flux  rope  eruptions  tend  to  show  stronger  writhe 
at  low  heights  than  ejective  eruptions  (CMEs).  This  argues  against  suggestions  that  the  writhing  facilitates  the  rise  of  the  rope  through 
the  overlying  field.  (4)  Erupting  filaments  which  are  S  shaped  already  before  the  eruption  and  keep  the  sign  of  their  axis  writhe  (which 
is  expected  if  field  of  one  chirality  dominates  the  source  volume  of  the  eruption),  must  reverse  their  S  shape  in  the  course  of  the  rise. 
Implications  for  the  occurrence  of  the  helical  kink  instability  in  such  events  are  discussed.  (5)  The  writhe  of  rising  loops  can  easily  be 
estimated  from  the  angle  of  rotation  about  the  direction  of  ascent,  once  the  apex  height  exceeds  the  footpoint  separation  significantly. 
Conclusions.  Writhe  can  straightforwardly  be  computed  for  numerical  data  and  can  often  be  estimated  from  observations.  It  is  useful 
in  interpreting  S  shaped  coronal  structures  and  in  constraining  models  of  eruptions. 

Key  words.  Magnetic  fields  -  Magnetohydrodynamics  (MHD)  -  Sun:  corona  -  Sun:  filaments,  prominences  -  Sun:  coronal  mass 
ejections  (CMEs) 


1.  Introduction 


Observations  of  the  solar  corona  display  a  variety  of  structures 
that  appear  S  shaped  when  viewed  in  projection  on  the  disk,  as 
for  example  filaments  with  curved  ends,  soft  X-ray  sigmoids, 
and  magnetic  loops  that  connect  different  active  regions.  The 
presence  of  an  S  shape  is  regarded  to  be  evidence  for  current- 
carrying  twisted  or  sheared  magnetic  fields  which  possess  mag¬ 
netic  helicity.  The  link  to  helicity  is  underlined  by  the  fact  that 
the  orientation  of  the  S  shows  a  hemispheric  preference,  as 
other  indicators  of  magnetic  chirality  (the  helicity  sign)  do  (e.g., 
iRust  &  Kumartl  19961:  IZirker  et  ai1ll997l).  Coronal  currents  store 
the  energy  required  to  power  eruptions  (lForbesll2000l).  Indeed, 
it  has  been  shown  that  active  regions  exhibiting  a  sigmoidal 
morphology  are  more  likely  to  erupt  than  non- sigmoidal  ones 
dCanfield  et  a1.l[T999h. 


Twisted  or  sheared  coronal  fields  carrying  nearly  force- 
free  volume  currents,  especially  magnetic  flux  ropes,  are 
central  in  models  of  filaments  and  prominences  (e.g., 
lAulanier  &  Demoulin  Il998l:  [Bobraetall  120081),  of  sig¬ 


moids  (e.g..  IRust  &  Kumar  Il996|:  iTitov  &  Demoulinl  Il999l: 
iLow  &  Berger!  120031:  iGibson  et  al.1  2004:  iKliem  et  alJ  120041). 
and  of  eruptions  (e.gjForbes  &  Isenberglll99ll:lAntiochos  et  al.l 


19991:  Amari  et  all  12000:  iFan  &  Gibsonl  120031:  iTorok  &  Klieml 
l  Yeates  &  Mackavl  20091).  The  amount  of  twist  or  shear  in 


2005 


them  can  be  quantified  by  the  magnetic  helicity.  For  a  magnetic 
flux  rope,  the  helicity  is  proportional  to  the  sum  of  its  twist  and 
writhe.  The  twist  measures  how  much  the  field  lines  wind  about 
the  magnetic  axis  of  the  rope,  whereas  the  writhe  quantifies  the 
helical  deformation  of  the  axis  itself. 


The  orientation  of  the  S  in  helical  structures  on  the  Sun  is 
highly  correlated  with  the  prevailing  sign  of  magnetic  helicity  in 
their  source  volume.  For  example,  IRust  &  Martin  (1994)  found 
a  one-to-one  correlation  between  the  chirality  of  sunspot  whirls 
and  the  orientation  (sinistral  vs.  dextral)  of  the  axial  field  of  fila¬ 
ments  that  spiral  into  their  periphery.  Here,  an  apparent,  or  true, 
clockwise  (counterclockwise)  rotation  of  the  sunspot,  associated 
with  positive  (negative)  helicity  of  its  field,  corresponds  to  for- 
ward  (reverse)  S  shape  of  the  filament  end.lPevtsov  et  al.l  (19971) 
found  that  forward  (reverse)  S- shaped  sigmoids  in  active  regions 
were  formed  in  predominantly  positive  (negative)  helicity  re¬ 
gions  for  «  90  %  of  the  cases  investigated.  In  the  remaining  cases 
the  helicity  sign  of  the  host  active  region  was  ambiguous,  which 
is  consistent  with  the  possibility  of  a  universal  relation  between 
the  orientation  of  the  S  in  such  structures  and  their  chirality.  In 
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Fig.  2.  The  writhe  of  a 
coil  depends  on  its  as¬ 
pect  ratio  (height/width) 
as  well  as  how  many 
turns  it  has.  Both  coils 
shown  have  3  turns,  but 
the  fat  coil  in  (a)  has  a 
writhe  W  =  2.68,  while 
the  thin  coil  in  (b)  has  a 
writhe  W  =  0.46. 


Fig.  1.  A  closed  flux  tube  as  seen  from  three  different  angles. 
The  number  of  apparent  crossings  one  sees  depends  on  the  view¬ 
ing  angle.  Crossings  can  be  called  left-handed  (negative)  or 
right-handed  (positive),  depending  on  the  relative  orientations  of 
the  upper  and  lower  field  lines.  The  diagrams  shown  in  (d)  give 
the  two  possible  cases  (up  to  rotation).  The  tube  shown  has  an 
average  crossing  number,  or  writhe ,  of  W  =  0.566.  Most  view¬ 
ing  angles  display  a  positive  crossing,  as  in  (b).  However,  some 
angles  do  not  show  any  crossings,  as  in  (a),  and  there  is  a  small 
set  of  angles  which  display  a  negative  crossing,  as  in  (c).  (The 
latter  picture  of  the  tube  is  obtained  by  viewing  it  from  the  top.) 


a  study  of  rotating  erupting  filaments.  1  Green  et  al.  d 2 0071)  found 
a  one-to-one  relationship  between  the  sign  of  dominant  helic- 
ity,  the  orientation  of  associated  sigmoids,  and  the  orientation  of 
the  rising  filament’s  developing  S  shape  (but  see  lMuglach  et  al.1 
1120091  for  occasional  exceptions  in  the  latter  association). 

In  addition  to  stable  structures,  erupting  filaments  and  promi¬ 
nences  can  also  exhibit  helical  deformations.  These  eruptions 
are  seen  as  rising  loops  in  the  corona,  which  in  most  cases  be¬ 
come  the  core  of  a  coronal  mass  ejection  (CME),  or  else  fall  back 
in  a  confined  eruption.  Their  change  of  shape  can  be  described 
as  a  writhing  of  the  filament  axis  out  of  the  plane  defined  by 
the  filament’s  foot  points  and  apex,  or  equivalently,  as  a  rota¬ 
tional  motion  of  the  filament  axis  in  the  upper  part  of  the  loop 
about  the  direction  of  ascent.  The  total  writhe  varies  widely  from 
event  to  event,  and  is  often  acquired  to  a  large  part  already  in  the 
low  and  middle  corona.  Some  very  clear  cases  of  such  writhing 
have  been  described  by  |Ji  etaD  020031).  iRomano  et  al]  (120031). 
IWilliams  et  al]  (120051).  andlZhou  et  al]  d2006l).  These  events  in¬ 
clude  a  wide  range  of  dynamical  behaviour,  from  a  confined 
eruption  to  the  fastest  CME  so  far  recorded.  The  writhing  in  the 
course  of  eruptions  has  attracted  interest  for  three  main  reasons. 

First,  it  indicates  that  the  erupting  field  has  the  magnetic 
structure  of  a  freely  moving  flux  rope,  line  tied  only  at  its  ends,  at 
the  onset  of  the  helical  deformation  (e.gjRus  I2OO3I: [Green  et  al .1 
120071). 

Second,  the  observed  deformations  correspond  exactly  to 
the  evolution  of  the  helical  kink  instability  of  the  current  chan¬ 
nel  in  the  core  of  a  flux  rope  (hereafter  KI)  and  have,  there¬ 
fore,  been  taken  as  strong  indication  of  this  instability’s  oc- 


currence  (Sakurai  19761;  Romano  et  al.  2003;  Rust  &  LaBonte 

2005 

:  Torok  &  Klierd  20051:  (Williams  et  al.1  2005b  Zhou  et  al. 

2006 

;  Gilbert  et  aLl(2007lk  The  KI  transforms  some  of  the  twist 

of  the  field  lines  about  the  magnetic  axis  of  the  rope  into  writhe 
of  the  axis.  In  doing  so,  it  lowers  the  magnetic  energy  by  reduc¬ 
ing  the  tension  in  the  twisted  field.  Since  this  is  an  ideal  MHD  in¬ 
stability  and  the  corona  is  a  nearly  perfectly  conducting  medium, 
the  conversion  of  twist  into  writhe  is  constrained  by  the  approxi¬ 
mate  conservation  of  the  magnetic  helicity  contained  in  the  rope 


(iBergerll  19841).  It  has  been  suggested  that  the  writhing  of  the  flux 
rope’s  upper  part  into  the  direction  of  the  overlying  field  is  en¬ 
ergetically  favourable  for  its  passage  through  the  overlying  field 
to  become  a  CME  dSturrock  et  al]  1200 It  iFanl  120051).  Thus,  the 
deformations  yield  hints  on  the  physical  processes  of  CME  initi¬ 
ation. 

Third,  the  writhing  is  the  major  factor  in  determining  the 
final  magnetic  orientation  of  the  CME  (apart  from  influences 
during  the  interplanetary  propagation),  which,  in  turn,  is  one  of 
the  critical  parameters  that  control  the  strength  of  the  interac¬ 
tion  if  the  CME  hits  the  terrestrial  magnetosphere — its  geoeffec¬ 
tiveness.  The  final  orientation  can  differ  largely  from  the  erupt¬ 
ing  structure’s  initial  orientation:  rotations  as  large  as  160°  have 
been  reported  (albeit  the  estimated  angles  are  very  approximate, 
seelDemoulin  2008;  Yurchvshvn  2008,  and  references  therein). 

Despite  the  common  occurrence  of  both  stable  and  dynamic 
helical  structures  on  the  Sun,  measurements  of  the  axis  writhe 
of  such  structures  have  not  yet  been  undertaken.  So  far,  only  es¬ 
timates  based  on  qualitative  considerations  have  been  reported 
(e.g.JVrsnak  et  al.lll993l:lRustll2003l).  This  is  mainly  due  to  the 
fact  that  the  observations  have  been  limited  to  a  2D  projection  of 
intrinsically  3D  structures  onto  the  plane  of  the  sky.  With  the  ad¬ 
vent  of  the  STEREO  mission,  true  3D  reconstructions  of  coronal 
structures  are  possible,  as  long  as  the  angular  separation  between 
the  twin  satellites  or  between  one  of  them  and  the  Earth  stays 
within  certain  intervals.  Even  for  numerical  simulations  which 
provide  the  3D  magnetic  field  of  evolving  unstable  flux  ropes, 
the  temporal  evolution  of  twist  and  writhe  has  only  extremely 
rarely  been  quantified  dLinton  et  al.ll  19981).  This  is  mainly  due 
to  the  fact  that  the  calculation  of  these  quantities  based  on  their 
general  definitions  is  a  relatively  complicated  task.  The  compu¬ 
tation  of  the  writhe  is  greatly  facilitated  by  jts  decomposition 
into  local  and  nonlocal  components  in  [Berger  &  Priori  (120061). 
Here  we  make  use  of  these  expressions. 

Our  aim  in  this  paper  is  to  describe  and  discuss  possible  ap¬ 
plications  of  writhe  measurements  for  both  stable  and  erupting 
objects  in  the  solar  corona.  To  this  end,  we  first  review,  in  Sect. [3 
the  concepts  of  twist  and  writhe  and  the  decomposition  of  the 
latter.  We  study  the  dependence  of  writhe  on  the  geometrical 
properties  of  helical  curves,  paying  particular  attention  to  their 
height  above  a  plane  and  to  the  presence  of  dips.  The  relation¬ 
ship  between  shape,  writhe,  and  helicity  of  S-shaped  flux  ropes 
is  then  discussed  (Sect.  [3]) .  In  Sect.  0  we  measure  the  conversion 
of  twist  into  writhe  in  several  numerical  simulations  of  unstable 
magnetic  flux  ropes  and  discuss  the  implications  of  the  results 
for  erupting  filaments  and  CMEs.  We  also  point  out  that  esti¬ 
mates  of  the  writhe  of  erupting  filaments  can  easily  be  obtained 
from  their  apex  rotation  as  soon  as  the  apex  height  significantly 
exceeds  the  footpoint  separation.  Section[5]  summarizes  the  find¬ 
ings. 
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Fig.  3.  A  loop  with  one  maximum,  seen  from  the  side  (a)  and 
the  top  (b).  The  nonlocal  writhe  is  proportional  to  the  difference 
in  angle  between  the  tangent  vector  at  the  top,  and  the  line  con¬ 
necting  the  footpoints.  This  angle  is  6h  -  #o  =  -7r/3  for  the  loop 
shown.  The  corresponding  nonlocal  writhe  is  Wnoniocai  =  1/3 
while  the  local  writhes  are  Wuocai  =  W2  local  -  -0.0573  for  a  to¬ 
tal  W  =  0.219.  (c):  The  writhe  for  the  family  of  curves  described 
by  Eqs.  GHI3.  The  horizontal  axis  gives  the  maximum  height 
h  of  the  curves,  measured  in  units  of  the  footpoint  separation. 
For  the  curves  shown,  the  nonlocal  writhe  Wnoniocai  =  -®/tt  is 
positive.  All  these  curves  display  a  reverse  S  shape  as  seen  from 
above. 


2.  Twist  and  Writhe 

2.1.  Definitions 

ICalugareanul  (1959)  introduced  a  quantity  called  writhe  to  mea¬ 
sure  how  much  a  closed  curve  coils  and  supercoils.  This  quan¬ 
tity  is  extensively  used  in  DNA  research  to  describe  the  physical 
configuration  of  the  DNA  molecule  dMaggioni  &  Riccal  12008) 
(in  addition  to  the  basic  double  helix  structure,  the  molecule 
must  be  highly  coiled  in  order  to  fit  into  a  microscopic  cell). 
There  are  several  equivalent  definitions  of  writhe  for  closed 
curves  (e.g.lAldinger  et  all  1 9951: iBerger  &  Priorll2006h.  Perhaps 
the  simplest  definition  to  understand  without  getting  lost  in 
mathematical  details  involves  counting  crossings  in  pictures  of  a 
curve  (see,  for  example,  Fig.  [I]).  If  we  draw  a  two-dimensional 
projection  of  the  curve,  it  crosses  itself  a  certain  number  of  times. 
Highly  coiled  or  tangled  curves  will  show  many  crossings,  while 
a  non-intersecting  curve  confined  to  a  plane  will  show  no  cross¬ 
ings  at  all  (unless  seen  exactly  edge-on).  Crossings  can  be  la¬ 
belled  positive  or  negative  depending  on  the  orientation  of  the 
sections  of  the  curve  above  and  below  the  crossing,  as  shown  in 
the  figure.  Thus  if  an  arrow  pointing  along  the  direction  of  the 
lower  curve  rotates  clockwise  to  match  the  direction  of  the  upper 
curve,  then  the  sign  is  positive.  Note  that  the  sign  of  a  crossing 


stays  the  same  if  we  reverse  the  direction  of  the  curve,  because 
both  segments  above  and  below  a  crossing  change  direction. 

Given  a  curve  and  a  viewing  angle,  we  can  count  the  num¬ 
ber  of  positive  crossings  N+  and  subtract  the  number  of  negative 
crossings  N- .  The  result  is  an  integer  -  however,  this  integer  de¬ 
pends  on  which  viewing  angle  we  are  using.  The  writhe  averages 
N+  -  N-  over  all  viewing  angles. 

As  a  consequence  the  writhe  of  a  circular  coil  depends  not 
only  on  how  many  coils  there  are,  but  on  whether  the  coils  are 
fat  or  thin,  i.e.  on  the  ratio  between  its  height  and  width  (Fig. [2]). 
For  thin  coils,  only  special  viewing  angles  (from  near  the  axis) 
will  see  crossings.  But  for  fat  coils  most  viewing  angles  display 
crossings. 

While  the  writhe  is  a  property  of  a  single  curve,  the  defini¬ 
tion  of  twist  involves  two  curves,  which  can  be  viewed  as  the 
edges  of  a  ribbon. ICalugareanul  C 1959)  derived  a  remarkable  for¬ 
mula  for  the  structure  of  a  ribbon.  Consider  two  closely  aligned 
curves,  like  the  two  sides  of  a  ribbon  or  the  two  strands  of  a 
DNA  molecule.  Assume  that  the  ribbon  or  molecule  closes  upon 
itself,  so  that  the  two  curves  do  not  have  endpoints.  These  two 
curves  link  each  other  by  some  (integer)  amount  L;  they  also 
twist  about  each  other  by  some  (usually  non-integer)  amount  T. 
Fet  the  writhe  W  be  the  writhe  of  one  of  the  two  curves  (alter¬ 
natively  a  central  curve  between  the  two).  Then  linking  number 
equals  twist  plus  writhe, 

L  =  T  +  W.  (1) 

For  our  purposes  we  need  two  modifications.  First,  we  are  con¬ 
cerned  with  magnetic  flux  tubes  rather  than  ribbons.  Also,  our 
tubes  are  not  closed:  they  have  ends  (footpoints)  on  a  boundary 
plane  or  planes.  For  these  two  reasons,  we  use  modified  expres¬ 
sions  for  writhe  as  given  in  IBerger  &  Prior  (  2006).  For  a  mag¬ 
netic  field,  we  can  replace  linking  number  L  with  magnetic  he- 
licity  //,  which  averages  the  linking  of  all  pairs  of  field  lines.  For 
a  straight  flux  tube  of  flux  F,  the  helicity  is  H  -TF2 .  (By  defi¬ 
nition,  if  field  lines  twist  around  the  axis  by  an  angle  of  O  =  2n, 
then  T  =  1 .)  The  twist  T  measures  how  much  of  the  helicity  is 
generated  by  parallel  electric  currents.  In  particular,  if  s  denotes 
arc  length  along  the  central  field  line  of  the  tube,  and  J\\  the  com¬ 
ponent  of  current  parallel  to  this  field  line,  then  T  -  f  ^ds,  with 

dT  _  hqJ\\ 
dv  4nB\\ 

The  writhe  for  a  curve  with  endpoints  (e.g.  the  axis  of  an  arbi¬ 
trary  shaped  flux  tube)  is  defined  by  the  magnetic  Calugareanu 
formula  (IBerger  &  Fieldlll984l:lMoffatt  &  Riccall  1 992l) 

H  =  (T  +  W)F2.  (3) 

A  current- free  tube  has  zero  twist,  by  Eq.  ([2]).  Thus  we  can  define 
the  writhe  of  a  curve  also  as  the  helicity  (divided  by  F2)  of  a  thin 
current- free  flux  tube  which  follows  the  curve. 

2.2.  Local  and  nonlocal  writhe 

iBerger  &  Prior  ( 2006)  give  expressions  for  the  writhe  of  curves 
stretching  between  two  parallel  planes,  as  well  as  loops  with 
both  endpoints  on  a  bottom  plane.  The  method  consists  of  sep¬ 
arating  the  curve  into  pieces  at  maxima  and  minima  in  height  z. 
Fet  us  suppose,  for  simplicity,  that  there  is  just  one  maximum, 
at  height  z  =  h.  We  first  ask  how  much  the  two  pieces,  or  legs, 
rotate  about  each  other  while  rising  from  the  bottom  to  the  top 
(see  Fig.  [3]).  Fet  6h  be  the  orientation  of  the  tangent  vector  at  the 
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Fig.  4.  The  writhe  for  a  family  of  dipped 
curves  with  height  function  given  by  a 
cubic  spline  (Eqs.  [TTHI41).  All  curves 
have  ) a  =  0.9.  (a):  the  height  function 
z(s)  with  maximum  height  h  -  0.6  for 
w  =  0.3, 0.4, 0.5, 0.6.  (b):  A  loop  with 
h  -  0.6,  w  =  0.4,  and  0  =  -n/3. 
The  projection  of  this  loop  onto  the  bot¬ 
tom  plane  is  the  same  as  in  Fig.  [3J),  but 
here  W  =  -0.03.  (c):  Writhe  calculated 
as  a  function  of  maximum  height  h  for 
fixed  0  =  -nil.  From  top  to  bottom, 
the  curves  have  w  =  0.3, 0.4, 0.5, 0.6 
(the  colors  correspond  to  panel  [a]),  (d): 
Writhe  calculated  as  a  function  of  max¬ 
imum  height  h  for  fixed  w  =  0.4.  From 
top  to  bottom  on  the  left,  the  curves  have 
0  =  — tt/8,  — tt/4,  — 3tt/8,  and  -n/2.  All 
these  curves  display  a  reverse  S  shape  as 
seen  from  above. 


top  of  the  loop  with  respect  to  the  x  axis.  Also  suppose  that  the 
line  from  the  positive  endpoint  in  the  bottom  plane  (z  =  0)  to  the 
negative  endpoint  has  an  orientation  0q-  Then  the  two  legs  of  the 
curve  rotate  about  each  other  by  6h  -  6q.  The  quantity 

W nonlocal  ~  (@h  ~  ^o)  (4) 
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(x(s),  y(s),  z(s))  for  0  <  v  <  1 .  We  choose  the  horizontal  coordi¬ 
nate  functions  separately  from  the  height  function.  A  family  of 
S  shapes  can  be  generated  using  the  form 


x(s) 

=  (s  -  1/2)  cos#(v); 

(7) 

y(s) 

=  (s  -  1/2)  sin^(v); 

(8) 

9(s) 

=  40j(1  -  s). 

(9) 

contributes  to  the  total  writhe.  It  is  called  the  nonlocal  writhe , 
because  for  most  of  their  lengths,  the  two  legs  are  far  away  from 
each  other.  Similar  formulae  can  be  derived  for  loops  with  more 
than  one  maximum  dBerger  &  Priorll2006i). 

In  addition,  each  leg  on  its  own  may  contribute  to  the  total 
writhe.  An  individual  leg  may  have  a  helical  shape,  for  example. 
The  individual  contributions  will  be  called  local  writhe.  Thus  for 
a  loop  with  legs  1  and  2,  we  have  the  decomposition 


W  —  Wi  local  T  W2 local  W nonlocal' 


(5) 


Let  Ti(z)  be  the  (unit)  tangent  vector  to  leg  1  at  height  z,  and  let 
T'  =  dTi/dz  be  its  derivative  with  respect  to  height.  Then  one 
finds 


1  (*  1 

Wl  local  -  2n  Jo  (1  +  |7Z|) 


(T,  xT'^dz, 


(6) 


with  a  similar  expression  for  W2  (seelBerger  &  Priodl 20061). 

If  only  the  S  shaped  projection  of  the  loop  onto  the  bottom 
plane  is  available,  as  it  is  often  the  case  in  solar  observations,  we 
can  still  infer  the  sign  of  the  local  writhe.  If  the  loop  bends  in 
a  clockwise  (counterclockwise)  manner  from  the  apex  to  either 
footpoint,  the  local  writhe  is  positive  (negative). 


2.3.  Height  dependence  of  writhe 

To  illustrate  how  the  writhe  behaves,  we  consider  a  family  of 
basic  loop  shapes.  Place  the  endpoints  of  the  loops  at  (x,y,z)  = 
(-1/2, 0,0)  and  (x,y,z)  =  (1/2, 0,0).  Thus  the  horizontal  dis¬ 
tance  between  the  endpoints  is  one;  equivalently,  all  lengths  will 
be  scaled  to  this  distance.  Let  the  axis  of  the  loop  be  a  curve 


The  number  0  determines  the  maximum  amount  of  rotation  of 
the  curve:  for  0  =  0  the  curve  remains  within  the  x  —  z  plane, 
while  for  0  =  ±n/ 2  the  top  of  the  curve  is  oriented  perpendicu¬ 
larly  to  the  line  between  the  endpoints. 

We  also  need  a  height  function  z(s).  Figure  [3^,b  employs  a 
parabolic  shape 

z(s)  =  4hs(l  -  s ),  (10) 

where  h  is  the  maximum  height,  and  0  =  — zr/3.  The  curve  ex¬ 
hibits  a  reversed  S  shape  when  seen  in  projection  on  the  bottom 
plane  (Fig.[3j3). 

For  curves  of  this  form,  local  and  nonlocal  writhe  give  con¬ 
tributions  of  opposite  sign.  We  can  ask  what  the  writhe  will  be 
as  a  function  of  height  h  and  maximum  rotation  0  (see  Fig.  [St). 
For  tall  curves,  the  nonlocal  writhe  dominates,  while  for  short 
curves  the  local  writhe  dominates.  Thus,  for  a  given  0,  the  sign 
of  writhe  depends  on  the  height  of  the  curve.  For  the  family  of 
curves  shown  in  Fig.  [3l  the  writhe  vanishes  for  all  curves  with 
h  ~  0.37.  Hence  the  writhe  can  be  zero,  even  if  the  curve  exhibits 
an  S  shape  when  seen  from  above. 

Changing  the  height  function  seems  to  make  little  difference, 
as  long  as  the  height  function  has  only  one  maximum.  For  exam¬ 
ple,  for  z(s)  =  h  sin(7rs)  we  find  almost  identical  values  for  the 
writhe;  the  crossover  height  where  writhe  vanishes  goes  down  to 
near  h  «  0.36.  A  quadratic  function  z(s)  =  h(  1  -  (2s  -  l)4)  gives 
similar  results,  with  crossovers  h  «  0.37-0.38. 

The  writhe  graphs  do  change  if  the  curve  has  more  than  one 
maximum.  Consider  a  loop  with  a  dip  at  its  central  part.  We  keep 
the  same  horizontal  coordinate  functions  as  defined  in  Eqs.  dTJ- 
0,  but  change  the  height  function  to  a  cubic  spline.  The  spline 
is  defined  by  the  conditions 

z(0)  =  z(l)  =  0; 


(11) 
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Fig.  5.  Heights  of  the  maxima  of  dipped  curves  at  which  the 
writhe  changes  sign,  for  different  dip  strengths  /t,  shown  as  a 
function  of  the  distances  between  the  maxima  of  the  height  func¬ 
tion  (see  text  for  details),  for  0  =  -7r/2.  Note  that  the  crossover 
heights  are  practically  independent  of  0  (Fig.  |4ji). 


z(0.5  -  w/2) 

=  z(0.5  +  w/2)  =  h; 

(12) 

z'( 0.5  -  w/2) 

=  z'(  0.5)  =  z\ 0.5  +  w/2)  =  0; 

(13) 

£(0.5) 

=  i uh; 

(14) 

in  addition,  second  derivatives  match  at  the  extrema.  The  param¬ 
eter  /i  gives  the  amount  of  dip,  e.g.  /i  =  0.9  gives  a  10  per  cent 
dip  (Fig.  0]),  and  w  is  the  distance  of  the  maxima.  Note  that  the 
horizontal  shape  has  not  changed;  the  dip  is  not  visible  if  the 
curve  is  seen  from  above  (Fig.0t>) . 

For  these  curves,  there  are  two  maxima  and  one  minimum. 
Let  the  tangent  vector  orientations  at  the  maxima  be  6\  max  and 
#2 max,  with  orientation  at  the  minimum  0m jn.  Then  we  replace 
Eq.  ®  with 

Wnonlocal  ~  (^lmax  T  ^2  max  —  ^min  —  $())•  (15) 

71 

For  such  curves,  the  writhe  often  has  the  same  sign  for  all  heights 
h  (see  the  bottom  two  curves  in  Fig.  0]:).  Note  that  forward 
(reverse)  S  shape  then  gives  positive  (negative)  writhe,  for  all 
heights.  For  some  dip  shapes  (e.g.,  the  top  two  curves  in  Fig. 0k) 
there  is  still  a  change  in  sign  of  writhe,  but  at  different  (and  usu¬ 
ally  larger)  heights  than  in  the  curves  without  a  dip.  For  exam¬ 
ple,  Fig.  0]i  shows  the  particular  case  w  =  0.4  for  four  different 
values  of  rotation  angle  0.  As  in  the  curves  without  dips,  the 
crossover  height  (here  near  h  =  0.9)  is  almost  completely  insen¬ 
sitive  to  the  rotation  angle. 

Figure  [5]  shows  that  the  range  of  w,  which  yields  a  sign 
change  of  writhe,  decreases  with  increasing  dip  strength  /i.  For 
w  >  0.45,  the  writhe  does  not  change  sign  for  all  /i  <  0.95. 
Hence,  in  order  for  a  sign  change  of  writhe  at  a  certain  height  to 
occur,  the  dip  must  be  relatively  small,  both  in  width  and  depth. 

3.  Sign  of  writhe  and  helicity  for  S-shaped  coronal 
flux  ropes 

The  height  dependence  of  writhe  for  loop- shaped  curves  ob¬ 
tained  in  Sect.  [2|  has  implications  for  the  relation  between  S 
shape,  sign  of  axis  writhe,  and  chirality  of  coronal  flux  ropes, 
which  we  discuss  in  the  following  two  sections. 


Since  a  projected  S  shape  of  a  curve  rooted  at  both  ends  in 
the  photosphere  is  an  immediate  signature  of  the  writhing  of  the 
curve  out  of  a  plane,  one  is  tempted  to  expect  a  unique  relation¬ 
ship  between  the  orientation  of  the  S  and  the  sign  of  the  curve’s 
total  writhe.  However,  as  shown  in  iBerger  &  Prior  (  2006)  and 
in  Sect.  12.31  this  relationship  is  not  unique.  The  ambiguity  can 
affect  the  relationship  between  the  orientation  of  S-shaped  flux 
ropes  and  their  chirality,  since  the  total  axis  writhe  contributes  to 
the  helicity.  This  leads  to  the  question  why  the  solar  observations 
indicate  such  a  strong  association  between  S  shape  and  chirality. 

Let  us  recast  the  two  striking  properties  of  curves  with  an 
S-shaped  projection,  which  are  apparent  from  the  results  in 
Sect.  12.31  First,  for  curves  without  a  dip,  the  sign  of  the  total 
writhe  will  flip  if  the  apex  height  h  is  changed  in  a  certain  range 
while  the  S  orientation  is  kept  (Fig.  [3]c).  Equivalently,  if  the  sign 
of  the  total  writhe  is  kept,  the  orientation  of  the  S  will  flip  (see 
Figs.lblandfTTlbelow).  The  flip  occurs  for  h  of  order  0.4  the  foot- 
point  distance,  which  is  relevant  for  solar  filaments.  Second,  the 
occurrence  or  disappearance  of  a  dip  (of  rather  moderate  depth 
[Fig.  01)  in  a  curve  of  sufficient  height  ( h  >  0.4)  can  flip  the  sign 
of  the  total  writhe  if  the  S  orientation  is  kept  (compare  Figs.  [3] 
and  0]). 

These  results  demonstrate  that  care  is  needed  when  deriving 
the  sign  of  the  writhe  from  the  observed  orientation  of  S-shaped 
structures  on  the  solar  disk.  It  must  be  ensured  that  the  height 
of  the  object  is  estimated  correctly,  and  the  possibility  that  the 
object’s  axis  has  a  dip  must  be  taken  into  account  if  h  >  0.4. 

First  consider  filaments.  The  height  of  stable  filaments  falls 
below  the  limit  h  «  0.4  in  most  cases.  The  presence  or  absence  of 
a  dip  is  then  irrelevant  for  the  relationship  between  the  S  orien¬ 
tation  and  the  sign  of  the  total  writhe.  For  erupting  filaments,  on 
the  other  hand,  it  is  often  possible  to  follow  their  evolution  into 
the  height  range  h  >  0.4,  especially  in  EUV  observations.  If  fila¬ 
ments  rise  to  such  and  greater  heights,  they  hardly  ever  display  a 
dip  in  their  top  part  (see  Sects.  l43l  and  14.41).  Thus,  in  practice,  S 
shape  and  writhe  can  be  unambiguously  related  to  each  other  if  a 
distinction  is  made  between  low-lying  (stable)  and  high-arching 
(erupting)  filaments  and  intermediate  cases  are  excluded. 

The  observed  correlation  between  the  magnetic  orientation 
(sinistral  vs.  dextral)  of  stable  active-region  filaments  and  the 
orientation  of  their  curved  ends  in  the  vicinity  of  sunspots  with 
whirls  dRust  &  Martini  1 994llZirker  et  al.ll  19971)  then  implies  that 
such  filaments  possess  positive  (negative)  writhe  if  the  orienta¬ 
tion  of  their  axial  field  is  sinistral  (dextral).  Thus,  the  sign  of 
writhe  is  identical  to  the  chirality  of  the  sunspot  field  near  the  end 
of  the  filament.  This  corroborates  the  expectation  that  such  fila¬ 
ments  are  threaded  by  field  emanating  from  the  sunspot,  which 
suggests  that  the  chirality  of  their  field  is  right  (left)  handed  for 
the  sinistral  (dextral)  orientation  of  the  axial  field.  Also,  it  allows 
for  the  possibility  that  writhe  helicity  contains  a  considerable 
fraction  of  the  helicity  of  the  field  that  threads  such  filaments. 

The  above  relationships  also  have  an  interesting  implication 
for  erupting  filaments  which  start  from  low  heights  ( h  <  0.4), 
possess  writhe  (S  shape)  already  before  they  rise  and  keep  the 
sign  of  writhe  during  their  evolution.  These  filaments  must  re¬ 
verse  their  S  orientation  in  the  course  of  the  rise.  We  will  con¬ 
sider  this  in  Sect.  14.31 

The  relation  between  S  shape  and  chirality  is  likely  to  be 
simpler  for  sigmoids.  These  sources  are  supposed  to  lie  in  sep- 
aratrix  surfaces  or  quasi- separatrix  layers  of  the  field  under¬ 
neath  flux  ropes  dTitov  &  Demoulinlll999h.  especially  if  they 
form  a  continuous  S  (  Gibson  et  al.  2004;  Green  &  Kliemll2009b. 
Consequently,  they  outline  low-lying  helical  field  lines  at  the  pe¬ 
riphery  of  a  flux  rope,  not  the  axis  of  a  rope.  The  writhe  of  these 
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Fig.  6.  (a,b):  Two  kink-unstable  flux  ropes,  displayed  by  isosurfaces  of  current  density  (from  Torok  et  al.  2004).  Both  simulations 
start  from  the  same  configuration  with  left-handed  twist  (i.e.  negative  helicity)  and  zero  axis  writhe,  but  with  oppositely  directed 
initial  perturbations.  The  top  panels  show  a  top  view  and  the  bottom  panels  show  a  side  view.  The  rope  in  (a)  kinks  downward  and 
develops  a  reverse  S  shape  when  seen  from  above.  The  rope  in  (b)  kinks  upward  and  develops  a  forward  S  shape  when  seen  from 
above,  (c):  Top  view  on  the  magnetic  axis  of  the  downward  kinking  rope,  shown  at  the  same  time  as  in  (a),  (d):  Same  as  (c)  for 
the  upward  kinking  rope.  Contour  plots  of  the  vertical  magnetic  field  in  the  bottom  plane  are  included  in  (c)  and  (d).  White  (black) 
contours  show  Bz  >  0  ( Bz  <  0). 


field  lines  has  the  same  sign  as  their  twist  (i.e.,  as  the  twist  helic¬ 
ity  of  the  rope),  and  it  has  a  unique  relationship  to  the  S  orienta¬ 
tion. 


4.  Writhe  of  erupting  flux  ropes  and  filaments 

In  this  section,  we  measure  the  writhe  of  the  magnetic  axis  of  a 
flux  rope  in  the  course  of  ideal  MHD  instabilities  in  numerical 
simulations,  and  we  discuss  the  implications  of  the  results  for 
filament  eruptions  and  CMEs.  In  flux  rope  geometry,  the  spine 
o£a  filament  follows  the  magnetic  axis  of  the  rope  closely  (e.g., 
lAulanier  &  Demoulinlll998l:TBobra  et  all!2008).  All  simulations 
presented  below  integrate  the  ideal  MHD  equations  and  use  the 
analytical  coronal  flux  rope  model  bvlTitov  &  Demoufinl  ([1999) 
as  initial  condition.  The  model  consists  of  a  line-tied,  arched, 
and  twisted  magnetic  flux  tube  embedded  in  an  arcade-like  po¬ 
tential  field  (seelTitov  &  Demoufinll  19991  for  details).  The  erup¬ 
tions  of  the  flux  rope  are  driven  either  by  the  helical  kink  insta¬ 
bility  (KI)  alone,  or  by  the  combined  action  of  the  KI  and  the 
torus  instability  dBatemanlll978l:lKfiem  &  Torokll2006l  hereafter 
TI).  Numerical  diffusion  permits  magnetic  reconnection  to  occur 
where  current  layers  steepen  in  response  to  the  development  of 
these  ideal  MHD  instabilities.  In  all  simulations,  we  measure  the 
writhe  of  the  rope’s  magnetic  axis  using  Eqs.  (l4l~[6landlT5l). 

4.1.  Kink-unstable  magnetic  flux  ropes 

Both  types  of  curves  discussed  in  Sect.  12.31  with  and  without 
a  dip  at  the  apex,  are  realised  in  the  3D  MHD  simulations  of 
kink-unstable  flux  ropes  bvlTorok  et  ali  (2004).  Figure  [6|  shows 
snapshots  of  the  current  channel  in  the  core  of  the  flux  rope  in 
the  course  of  the  instability  for  two  of  these  simulations.  Both 
simulations  start  from  the  same  flux  rope  configuration  with  a 
left-handed  average  twist  angle  of  O  =  2 nT  -  -4.9  n  and  van¬ 
ishing  writhe  of  the  axis  (whose  projection  on  the  bottom  plane 
is  a  straight  line).  In  the  following  we  will  use  the  twist  angle  in 
quantitative  statements.  Its  end-to-end  value  is  calculated  as  an 
average  over  the  cross-section  of  the  current  channel  in  the  core 
of  the  flux  rope  (iTorok  et  al.ll2004l). 

The  downward  kinking  rope  shown  in  Fig.  [6^  has  a  dip  in 
its  middle  part  and  displays  a  reverse  S  shape  when  seen  from 
above.  Using  the  parameters  defined  in  Sect.  12.31  its  axis  can  be 


described  by  h  -  0.26,  w  =  0.57,  p  =  0.60,  and  0  =  -7r/4.4. 
The  rope  axis  runs  relatively  flat,  has  a  strong  dip,  and  the  two 
maxima  are  relatively  far  away  from  each  other.  From  our  con¬ 
siderations  in  Sect.  12  31  we  therefore  expect  the  writhe  to  be  neg¬ 
ative.  The  negative  sign  of  the  writhe  follows  also,  of  course, 
from  the  conservation  of  magnetic  helicity  in  the  course  of  the 
instability.  The  axis  writhe  is  W  =  -0.26,  with  Wiocai  =  0.23, 
Wnoniocal  =  -0.49  (the  latter  being  made  up  of  contributions 
-0.23  at  the  minimum  and  -0.13  at  each  maximum). 

The  axis  of  the  upward  kinking  rope  develops  a  forward  S 
shape  when  seen  from  above  (Fig.  (6t>).  It  has  only  one  max¬ 
imum,  at  h  =  0.53.  This  height  is  larger  than  the  crossover 
height  for  curves  with  one  maximum  (Sect.  12.31).  hence  we  ex¬ 
pect  the  writhe  to  be  negative.  Again,  this  follows  from  helic¬ 
ity  conservation.  We  find  W  =  -0.30,  with  Wiocai  =  -0.01  and 
Wnoniocal  ~  — 0.29  (0  —  7T / 3 .5).  Note  that  Wiocai  tffld  Wnoniocai 
have  the  same  sign,  which  is  never  the  case  for  the  family  of 
S-shaped  curves  with  one  maximum  discussed  in  Sect.  12.31  This 
discrepancy  is  due  to  the  fact  that  the  flux  rope  axis  in  the  simula¬ 
tion  is  not  perfectly  S-shaped  when  seen  from  above  (Fig.[6t,d). 
Rather,  close  to  the  footpoints,  the  axis  bends  in  the  direction 
opposite  to  the  overall  orientation  of  the  S,  contributing  nega¬ 
tive  local  writhe.  This  additional  bending  is  an  indication  for  the 
occurrence  of  KI  eigenmodes  with  axial  wavenumbers  k  >  1 
(see,  e.g..  iLinton  et  al.lll998l).  which  seems  plausible  given  the 
relatively  large  twist  used  in  the  simulations.  As  a  consequence, 
the  local  writhe  changes  sign  as  one  follows  the  axis  from  the 
apex  to  either  footpoint,  which  is  not  the  case  for  the  perfectly 
S-shaped  curves  discussed  in  Sect.  12.31 

Note  that  the  writhe  of  the  two  kinked  ropes  is  very  simi¬ 
lar.  Thus,  the  two  simulations  illustrate  that  flux  ropes  of  equal 
writhe  can  have  opposite  orientations  of  their  projected  S-shape, 
depending  on  the  apex  height. 

Let  us  now  consider  the  amount  of  twist  that  is  converted 
into  writhe  in  the  course  of  the  KI  for  the  cases  shown  in  Fig.  [6] 
Since  the  KI  is  an  ideal  MHD  instability,  both  the  axial  flux  and 
the  helicity  of  the  flux  rope  must  be  conserved,  i.e.,  twist  must 
be  converted  into  an  equal  amount  of  writhe  of  the  same  sign. 
However,  we  must  take  care  when  calculating  the  twist  con¬ 
verted  within  the  entire  rope,  since  flux  surfaces  far  away  from 
the  axis  might  bend  in  a  different  way  than  the  axis  (for  example 
if  the  flux  rope  expands  significantly).  In  the  simulations  shown 
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Fig.  7.  Confined  filament  eruption  on  2002  May  27  observed 
by  TRACE  in  the  195  A  band,  and  magnetic  field  lines  outlining 
the  core  of  the  kink-unstable  flux  rope  in  the  confined  eruption 
simulation  by  Torok  &  Kliem  (2005).  The  flux  rope  axis  in  the 
simulation  has  a  writhe  of  W  =  0.67  at  the  state  shown  (t  =  37 r a 
[Alfven  times],  apex  height  =1.1  times  the  footpoint  distance; 
see  also  Fig.  [9^).  Note  that  the  flux  rope  twist  in  the  simulation 
was  chosen  right-handed,  in  order  to  account  for  the  morphology 
of  the  kinked  filament. 


in  Fig.[6l  this  is  not  the  case,  so  we  can  estimate  the  twist  con¬ 
verted  within  the  entire  rope  (not  just  in  the  vicinity  of  its  axis) 
from  the  axis  writhe.  The  writhe  W  =  -0.26  ( W  =  -0.30)  for 
the  downward  (upward)  kinking  rope  corresponds  to  a  converted 
twist  of  O  =  -0.52  n  (®  =  -0.60  n).  In  order  to  check  the  reli¬ 
ability  of  this  estimate,  we  independently  measure  the  flux  rope 
twist  for  the  upward  kinking  case  at  the  stage  of  the  evolution 
shown  in  Fig.  [6)3,  using  Eqs.  (12-13)  in  [Berger  &  Priori  (120061). 
We  find  that  the  converted  twist  measured  that  way  and  the  value 
estimated  from  the  axis  writhe  differ  by  less  than  1%,  as  long  as 
we  exclude  flux  surfaces  very  close  to  the  rope  surface,  where 
the  twist  is  strongly  nonuniform  (see  Fig.  2  inlTorok  et  all 20041). 
This  shows  that  the  twist  converted  in  the  course  of  the  KI  can 
be  reliably  estimated  from  the  axis  writhe  if  the  inhomogeneity 
of  the  radial  twist  profile  and  its  changes  remain  modest.  It  also 
shows  that  the  helicity  is  conserved  to  a  high  degree  of  accuracy 
in  the  simulation. 

Interestingly,  the  amount  of  converted  twist  of  «  (0.5  -0.6)  n 
at  the  stage  of  the  instabilities  shown  in  Fig.  [6]  is  smaller  (only 
«  10  %  of  the  initial  twist)  than  one  might  intuitively  expect 
from  the  relatively  strong  deformation  of  the  flux  ropes.  It  is 
often  assumed  that  the  KI  converts  a  twist  of  ±2  n,  which  proba¬ 
bly  arises  from  the  fact  that  helically  deformed  structures  on  the 
Sun  typically  exhibit  one  helical  turn,  and  are  therefore  often  be¬ 
lieved  to  have  a  writhe  of  ±1  fe.g.lRustll2003L  see  also  Sect.  4.2). 
However,  as  Fig.  [2]  shows,  the  number  of  helical  turns  of  a  curve 
or  flux  rope  can  be  very  different  from  its  writhe,  so  that  estima¬ 
tions  of  writhe  (and  of  twist  converted  in  the  course  of  the  KI) 
cannot  be  made  from  the  observed  number  of  turns  alone.  In  the 
following  subsection,  we  will  investigate  the  conversion  of  twist 
into  writhe  in  the  course  of  flux  rope  instabilities  in  more  detail. 

4.2.  Writhe  in  confined  and  ejective  eruptions 

[Torok  &  Klieml  (2005)  presented  two  simulations  of  erupting 
kink-unstable  flux  ropes  with  identical  initial  average  flux  rope 
twist,  I® |  ~  57T,  and  nearly  identical  geometrical  rope  parame¬ 
ters,  but  different  ambient  potential  fields.  In  the  first  simulation 
(hereafter  TK1;  Fig.  [7]),  the  potential  field  decays  only  slowly 
with  height  above  the  initial  flux  rope  position.  The  rope  ascends 
driven  by  the  KI  (similar  to  the  case  shown  in  Fig.  [SJd),  but  its 
rise  terminates  at  an  apex  height  of  «  3.5  times  the  initial  height 


Fig.  8.  (a):  Filament  eruption  on  2003  February  18  observed 
by  the  EUV  Imaging  Telescope  (EIT)  onboard  the  Solar  and 
Heliospheric  Observatory  in  the  195  A  band  at  02:12  UT.  (b): 
The  same  filament  as  in  (a),  observed  as  CME  core  by  HAO’s 
Mauna  Loa  Solar  Observatory  Mk4  white-light  coronagraph  at 
02:34  UT.  The  filament  as  seen  by  EIT  does  not  show  indica¬ 
tions  for  the  occurrence  of  a  significant  rotation  about  its  rise 
direction,  whereas  the  CME  core  does  (note  the  crossings  of 
its  legs),  (c):  Two  snapshots  of  the  CME  simulation  in  Torok 
&  Kliem  (2005)  at  the  same  viewing  angle  but  different  spatial 
scale,  taken  at  t  -  22r a  (left)  and  t  -  43ta  (right).  The  apex 
height  of  the  flux  rope  axis  is,  respectively,  0.88  and  3.8  of  the 
footpoint  distance  (which  is  «  3.3  times  the  initial  apex  height). 
The  total  (nonlocal)  writhe  of  the  axis  is,  respectively,  -0.254 
(-0.212)  and  -0.623  (-0.518). 


ho,  in  very  good  quantitative  agreement  with  a  confined  filament 
eruption  observed  on  2002  May  27  fseelJi  et  al.ill2003l). 

In  the  second  simulation  (hereafter  TK2;  Fig.  [8]),  the  poten¬ 
tial  field  decreases  much  faster  with  height.  As  a  result,  the  TI 
sets  in  once  the  KI  has  lifted  the  rope  to  a  height  where  the  po¬ 
tential  field  drops  off  sufficiently  steeply.  The  rope  is  then  addi¬ 
tionally  accelerated  by  the  TI  and  eventually  ejected.  The  flux 
rope  rise  characteristics  could  be  scaled  to  closely  match  the  ac- 
celeration  profile  of  a  CME  on  2001  May  15  (seelMaricic  et  al.l 
12004. 

Here  we  measure  the  evolution  of  the  writhe  of  the  flux  rope 
axis  in  these  simulations.  As  in  the  simulations  described  in 
Sect.  14.11  the  writhe  vanishes  initially.  Figured  shows  the  evolu¬ 
tion  of  W ,  of  W nonlocal,  and  of  the  apex  rotation.  The  rotation  is 
proportional  to  Wnoniocau  since  the  axis  has  only  one  maximum 
at  all  times.  The  writhe  first  grows  exponentially,  followed  by 
a  transition  to  a  saturation  phase.  The  second  strong  increase  in 
TK1  at  t  «  35ta  is  caused  by  an  additional  deformation  of  the 
flux  rope  axis  due  to  the  onset  of  magnetic  reconnection  with 
the  overlying  field  around  this  time  (seelTorok  &  Klieml 120051). 
Wnoniocal  reaches  «  0.5  in  TK2,  corresponding  to  a  rotation  of 
«  90°.  A  very  similar  value  is  found  for  the  rope  in  TK1  before 
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Fig.  9.  (a):  Absolute  values  of  total  writhe  (dashed)  and  nonlo¬ 
cal  writhe  (solid)  as  a  function  of  time  for  the  confined  (blue) 
and  ejective  (red)  flux  rope  eruption  simulations  in  Torok  and 
Kliem  (2005).  (b):  Absolute  values  of  flux  rope  axis  apex  ro¬ 
tation  in  degrees  for  the  two  simulations,  as  a  function  of  axis 
apex  height  normalized  to  its  initial  value.  Note  that  the  rota¬ 
tion  is  proportional  to  the  nonlocal  writhe  (a  nonlocal  writhe  of 
±0.5  corresponds  to  a  rotation  of  +90°).  We  plot  absolute  val¬ 
ues  here  since  different  signs  (or  handedness)  of  the  initial  flux 
rope  twist  were  used  in  the  simulations.  Time  is  given  in  Alfven 
times  based  on  the  initial  apex  height  of  the  flux  rope  axis  and 
the  initial  Alfven  velocity  at  this  point,  ta  =  ho/VAO- 


it  starts  to  reconnect  with  the  overlying  field.  As  expected  from 
Sect.  12.31  the  nonlocal  writhe  clearly  dominates  the  local  writhe 
once  the  flux  rope  has  sufficiently  risen.  In  both  simulations,  the 
rope  develops  a  clear  S  shape  when  viewed  in  projection  on  the 
bottom  plane. 

The  temporal  evolution  of  the  writhe  is  very  similar  in  TK1 
and  TK2,  since  the  growth  rate  of  the  KI,  largely  set  by  the  initial 
twist,  is  nearly  the  same  (Fig.  [9 k).  The  evolution  as  a  function 
of  height,  however,  is  quite  different  (Fig.  03).  The  additional 
acceleration  of  the  rope  by  the  TI  spreads  the  rotation  of  the 
apex  over  a  larger  height  range  (compare  also  Figs.  1  and  4  in 
[Torok  &  KfienJl2005h.  The  TI,  which  is  a  form  of  the  lateral 
kink  instability,  primarily  expands  the  unstable  flux  loop,  while 
the  KI  primarily  produces  a  helical  shape. 

Given  the  nearly  identical  choice  of  the  initial  flux  rope  pa¬ 
rameters  in  the  two  simulations,  one  can  conclude  from  Fig.  [9)3 
that  a  stronger  field  immediately  above  the  initial  rope  does  not 
only  resist  the  evolution  into  a  CME  more  efficiently,  but  also 
produces  a  more  pronounced  writhing  at  low  heights.  This  as¬ 
sociation  is_opposite  to  the  suggestion  in ISturrock  et  al.l  (2001) 
and  [Fan  (  2005)  that  the  writhing  facilitates  the  rupture  of  an  un¬ 
stable  flux  rope  through  the  overlying  field  (which  is  of  course 
strongest  at  low  heights),  and  it  underlines  the  importance  of 
magnetic  reconnection  below  the  rope  in  permitting  the  erupting 
flux  rope  to  pass  through  the  overlying  field  (e.gjLin  &  Forbesl 
120001: 1 V  r  snakH2008h . 

Since  the  writhing  in  CMEs  tends  to  be  distributed  over  a 
large  height  range  (Fig.  03),  much  of  it  may  escape  detection. 
Typically,  the  writhing  is  apparent  from  the  apex  rotation  of  an 
associated  erupting  filament  or  prominence  observed  in  H a  or 
in  the  EUV,  and  this  is  usually  limited  to  the  low  and  middle 
corona.  Therefore,  a  significant  rotation  (and  the  KI)  may  occur 
in  a  larger  fraction  of  CMEs  than  usually  thought. 

An  illustrative  example,  the  filament  eruption  and  CME  on 
2003  February  18  (lFanll2005l:  iRust  &  LaBontell2005l).  is  shown 
in  Fig.0  The  EIT  data  do  not  yield  indications  of  a  significant 
rotation,  but  the  Mk4  coronagraph  reveals  that  the  legs  of  the 
CME  core  cross  to  form  an  “inverse  y”  shape  when  the  core 
reaches  a  height  of  «  1 RQ  above  the  solar  surface.  Such  a  shape 


develops  due  to  writhing  and  is  commonly  regarded  to  be  evi¬ 
dence  of  the  KI  (e.g. .[Gilbert  et  al.ll2007l). 

We  note  that  the  inverse  y  shape  has  been  associated  for  the 
two  events  shown  in  Figs.  [7]  and  [8]  with  a  writhe  of  +1  and  —  1, 
respectively  dRustl  120031:  IRust  &  LaBontell2005l).  However,  the 
values  obtained  from  the  simulations  indicate  that  the  writhe  of 
such  structures  is  rather  some  number  between  0  and  1 .  This  ap¬ 
pears  plausible  if  one  recalls  that  the  writhe  can  be  expressed 
as  average  crossing  number  (see  Sect.  12.  Hand  Fig.Q]).  Consider, 
for  example,  the  filament  shown  in  Fig.  0  from  the  particular 
viewing  angle  of  the  observation  we  see  one  (positive)  crossing 
of  the  axis.  Other  viewing  angles  will  display  either  one  cross¬ 
ing  or  no  crossing,  so  that  the  average  crossing  number,  i.e.  the 
writhe,  must  be  less  than  unity. 

4.3.  Changing  S  shape 

The  simulations  presented  in  Sects.  I47TI  and  14.21  all  start  from  a 
flux  rope  which  is  straight  in  projection  on  the  bottom  plane. 
Filaments,  however,  often  exhibit  an  S  shape  already  in  their 
equilibrium  state.  In  general,  the  spine  of  such  filaments  is  a 
curve  of  nonvanishing  writhe.  In  considering  their  eruption,  we 
assume  that  the  magnetic  field  has  formed  a  flux  rope  topology 
before  or  in  the  early  stages  of  the  eruption.  If  the  filament  erupts 
from  a  low  height  ( h  <  0.4),  keeping  the  sign  of  writhe  (which 
is  expected  if  one  chirality  dominates  the  field  in  the  source  vol¬ 
ume),  then,  according  to  Sects.  0 and 0  a  change  of  the  S  orien¬ 
tation  must  occur. 

Figure  ffOli-b  shows  a  filament  eruption  in  which  an  initial 
strong  reverse-S  shaped  bending  is  completely  straightened  out. 
A  change  to  a  forward  S  may  have  followed,  suggested  by  the 
rapid  rotation  of  the  filament.  However,  the  eruption  occurred 
about  700  arcsec  from  disc  center,  so  that  the  expected  transition 
to  a  forward  S  could  not  be  witnessed.  Figure  lT0b-e  shows  an 
erupting  filament  that  formed  an  incomplete  S.  As  expected,  it 
changed  its  J  shape  from  reverse  to  forward  while  rising.  See 
iRomano  et  al.l  (120051).  IRust  &  LaBontel  (120051).  andlGreen  et  al.l 
(120071)  for  descriptions  and  discussions  of  the  apex  rotation  and 
the  magnetic  environment  in  these  events. 

Observationally,  it  is  difficult  to  confirm  the  dynamic  tran¬ 
sition  from  one  complete  S  shape  to  the  other  (e.gjGreen  et  al.l 
1 20071).  In  order  to  demonstrate  the  transition,  we  perform  a  simu¬ 
lation  similar  to  the  ones  described  in  Sect.  14.21  but  now  starting 
from  an  initially  S- shaped  flux  rope.  To  this  end,  we  construct  a 
stable  numerical  equilibrium  of  a  modified  version  of  the  Titov 
&  Demoulin  (T&D)  model,  which  contains  such  a  rope.  The 
rope  has  left-handed  average  twist  of  «  4.7  n,  similar  to  the  ini¬ 
tial  twist  in  the  simulations  described  above.  Note  that  the  helical 
shape  of  the  constructed  flux  rope  permits  a  larger  twist  in  stable 
equilibrium  than  the  original,  toroidal  T&D  rope  of  identical  as¬ 
pect  ratio  in  rather  similar  ambient  field.  We  trigger  the  eruption 
of  the  flux  rope  by  imposing  slow,  quasi-static  converging  flows 
towards  the  polarity  inversion  line  in  the  bottom  boundary.  The 
construction  of  the  modified  equilibrium  and  the  details  of  the 
simulation  will  be  described  in  a  separate  publication.  Here  we 
are  merely  interested  in  the  evolution  of  writhe  and  S  shape  of 
the  erupting  flux  rope. 

Initially,  the  flux  rope  has  a  dip  in  its  middle  part  and  dis¬ 
plays  a  clear  reverse  S  shape  when  seen  from  above  (Tig.fTTh). 
Using  the  parameters  discussed  in  Sect.12.31  its  axis  can  be  char¬ 
acterised  by  0  =  — tt/3.5,  h  =  0.17,  p  =  0.91,  and  w  =  0.43. 
Inspecting  Fig.  ®  we  expect  an  axis  writhe  of  about  -0.2  for 
these  values.  Measuring  the  axis  writhe,  we  find  W  =  -0.07. 
As  in  the  above  simulations,  the  writhe  is  smaller  than  expected 
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Fig.  10.  Top  row:  Filament  eruption  on  2001  June  15  observed 
by  TRACE  in  195  A.  (a):  At  09:57  UT,  just  before  the  eruption. 
The  filament  exhibits  a  clear  inverse  S  shape  in  its  upper  part, 
(b):  At  10:06,  during  the  eruption.  The  shape  of  the  filament 
has  straightened  out.  Flare  ribbons  have  formed.  Bottom  row, 
(c-e):  Filament  eruption  on  2000  June  6  observed  by  TRACE  in 
171 A  at  13:29:27, 13:31:34,  and  13:32:20  UT,  respectively.  The 
filament  (outlined  by  yellow  dots)  changes  its  orientation  from  a 
reverse  J  to  a  forward  J. 


Fig.  11.  Perspective  and  top  view  on  magnetic  field  lines  out¬ 
lining  the  flux  rope  core  in  the  simulation  described  in  Sect.  14.31 
(a):  after  the  initial  relaxation  (time  reset  to  t  -  0).  The  two 
maxima  of  the  flux  rope  axis  have  a  height  of  0.17  of  the  foot- 
point  separation  and  0  =  -n/3.5.  (b)  and  (c):  in  the  course  of 
the  eruption,  at  t  -  55ta  and  64ta,  respectively.  The  writhe  is 
W  =  -0.07,  0.01,  and  -0.15  (from  left  to  right). 


from  the  curves  discussed  in  Sect.  12.31  since  also  here  the  axis 
bends  opposite  to  the  overall  S  close  to  the  flux  rope  footpoints. 
However,  only  the  sign  of  the  writhe  is  important  for  the  present 
discussion. 

As  the  converging  flows  are  applied,  the  flux  rope  first  rises 
slowly  in  response  to  the  progressive  weakening  of  the  tension 
of  the  overlying  potential  field.  After  t  ~  45ta  it  accelerates 
rapidly,  most  likely  due  to  the  onset  of  the  KI  or  TI  (we  refrain 
here  from  a  detailed  investigation  of  the  acceleration  mecha¬ 
nism).  As  the  rope  rises,  its  shape  first  straigthens  out,  similar 
to  the  filament  eruptions  shown  in  Fig.  [TOj  and  the  axis  writhe 
goes  to  zero.  At  the  state  shown  in  Fig.  Eb  (t  =  55ta),  the  axis 


has  developed  a  single  maximum  (i.e.,  the  dip  has  disappeared), 
its  height  is  0.27,  and  its  writhe  is  W  =  0.01.  Subsequently,  the 
rope  rapidly  develops  a  forward  S  shape,  and  the  writhe  becomes 
negative  again.  At  the  state  shown  in  Fig.  Et  ( t  =  64ta),  one 
finds  h  =  0.69  and  W  =  -0.15.  Except  for  the  short  period  when 
the  writhe  is  close  to  zero  (as  is  expected  since  the  axis  becomes 
almost  straight  in  projection  in  the  reversal),  it  stays  always  neg¬ 
ative  during  the  rise  of  the  rope.  Therefore,  a  transition  from  a 
reverse  to  a  forward  S  shape  must  occur. 

A  qualitatively  similar  transition  occurred  in  the  simulation 
of  a  kink-unstable  erupting  flux  rope  by  iFanl  (  2005:  see  her 
Fig.  6a-b),  although  the  flux  rope  equilibrium,  the  external  po¬ 
tential  field,  and  the  driving  in  the  photospheric  boundary  all 
differed  from  our  simulation.  Note  that  Fan  displays  only  field 
lines  below  the  magnetic  axis  of  the  flux  rope,  but  these  extend 
up  to  the  vicinity  of  the  axis,  so  that  the  actual  rope  axis  must 
behave  in  the  same  manner.  Also  note  that  the  straightening  of 
the  S  occurs  in  the  slow-rise  phase  for  both  simulations.  These 
results  have  interesting  consequences  for  the  occurrence  of  the 
KI  in  filament  eruptions. 

First,  it  is  important  to  note  that  both  systems  are  constructed 
such  that  only  one  sign  of  helicity  is  present  at  the  onset  of  the 
rise.  Thus,  the  initial  reduction  of  the  writhe  (from  a  negative 
value  to  nearly  zero)  cannot  be  driven  by  a  conversion  of  flux 
rope  twist  into  writhe  under  conservation  of  helicity,  since  the 
twist  has  the  same  sign  as  the  initial  writhe.  In  other  words,  if  an 
equilibrium  flux  rope  possesses  writhe  (i.e.,  an  S  shape),  its  ini¬ 
tial  slow  rise,  which  must  reduce  the  S  shape  to  a  nearly  straight 
rope,  cannot  be  caused  by  a  helical  kink  instability,  since  this 
instability  transforms  twist  into  writhe.  A  different  process  is  re¬ 
quired,  which  actually  transforms  writhe  into  twist.  The  details 
of  this  process  will  likely  require  further  study,  but  it  is  obvious 
that  the  converging  flows  in  the  simulation  of  Fig.  \lT\  and  the 
continuing  emergence_of  a  twisted  flux  rope  through  the  bottom 
boundary  inlFanl  (120051)  were  the  drivers  in  the  two  simulations. 

Second,  the  ability  of  the  S-shaped  flux  rope  to  remain  in 
stable  equilibrium  for  a  higher  twist  than  a  straight  flux  rope 
could  accommodate,  and  the  writhe  conversion  into  twist  in  the 
slow-rise  phase,  both  support  the  occurrence  of  the  KI  after  the 
flux  rope  has  become  straight.  The  two  simulations  demonstrate 
a  scenario  that  enables  the  KI  in  a  system  which  is  driven  signifi¬ 
cantly  beyond  the  threshold  of  the  instability  by  the  straightening 
in  the  slow-rise  phase,  permitting  a  significant  growth  rate. 

4.4.  Writhe  estimate  for  rotating  eruptive  filaments 

The  dominance  of  the  nonlocal  writhe  over  the  local  writhe  for 
tall  curves  discussed  in  Sect.  [23] provides  the  possibility  of  esti¬ 
mating  the  writhe  of  rotating  filaments  that  erupt  near  Sun  cen¬ 
ter  even  if  no  3D  observations  are  available.  When  the  filament 
has  reached  a  sufficient  height,  the  writhe  can  be  approximated 
by  the  nonlocal  writhe  (see  Fig.  [9^).  The  latter  can  easily  be 
obtained  from  Eq.  @  by  inserting  the  observed  rotation  angle 
(Fig.  [3})).  For  this  estimate,  one  has  to  verify  that  the  erupting 
filament  does  not  possess  a  dip  in  its  upper  part  in  the  relevant 
height  range. 

Limb  observations  of  erupting  filaments  and  prominences  do 
not  indicate  the  presence  of  such  a  dip.  We  can  ask  under  which 
circumstances  a  central  dip  in  an  erupting  flux  rope  might  occur. 
In  the  simulation  described  in  Sect.  14.31  the  flux  rope  has  a  dip 
prior  to  eruption.  As  the  rope  rises,  the  dip  disappears  early  in 
the  evolution,  long  before  the  rope  reaches  a  height  where  the 
nonlocal  writhe  dominates  the  local  writhe  (Tig.fTTb).  A  dip  in 
the  central  part  of  a  flux  rope  will  naturally  develop  if  the  rope 
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Fig.  12.  Snapshots  of  two  simulations  of  kink-unstable  flux  ropes  with  an  initial  average  twist  of  O  =  —  1  l7r.  Magnetic  field  lines 
outlining  the  flux  rope  core  as  well  as  contours  of  the  normal  magnetic  field  in  the  bottom  plane  are  shown,  (a):  Initial  configuration, 
(b):  Upward  kinking  flux  rope,  (c):  Downward  kinking  flux  rope. 


is  kink-unstable  and  is  perturbed  such  that  its  apex  moves  down¬ 
ward,  as  in  the  simulation  shown  in  Fig.  [6^.  However,  in  this 
simulation  the  rope  does  not  erupt,  rather  its  central  part  is  con¬ 
tinuously  moving  downward  while  the  remaining  parts  of  the 
rope  do  not  ascend  significantly. 

It  is  well  known  that  a  kink-unstable  flux  rope  can  de- 
velop  more  than  one  helical  turn  if  it  is  sufficiently  twisted  (e.g. 
iLinton  et~aflll998l).  For  a  flux  rope  with  two  turns,  a  dip  might 
form  between  the  two  helical  turns.  To  check  this,  we  perform  a 
simulation  similar  to  the  ones  described  in  Sects. 14.  lland!4.2[  but 
now  we  consider  a  very  large  initial  flux  rope  twist  of  O  =  -lln 
(Fig.fT2k).  As  expected,  the  rope  undergoes  a  strong  helical  de¬ 
formation.  However,  no  dip  forms  in  its  central  part.  Rather,  two 
dips  develop  along  the  legs  of  the  rope,  while  its  central  part 
ascends  (Fig.fl2b).  For  such  a  configuration,  care  must  be  taken 
when  estimating  the  writhe  from  the  apex  rotation,  since  the  axis 
rotation  at  the  other  extrema  contributes  to  the  nonlocal  writhe 
as  well.  However,  such  astrong  writhing  is  only  very  rarely  ob¬ 
served  (see,  e.gjRomano  et  al 120031  for  a  filament  eruption  that 
corresponds  quite  well  to  the  simulation  in  Fig.fl2b). 

When  the  simulation  is  repeated  with  a  small,  downward  di¬ 
rected  initial  velocity  perturbation,  a  remarkably  different  shape 
is  obtained.  Now  the  rope  develops  a  large  dip  in  its  central  part. 
As  in  the  simulation  shown  in  Fig.  [6^,  the  central  section  moves 
downward.  However,  in  contrast  to  that  simulation,  the  remain¬ 
ing  parts  of  the  rope  now  rise  significantly  fFig.[l2b).  Still,  they 
do  not  succeed  in  pulling  the  central  part  upward,  that  part  rather 
keeps  moving  downward  until  it  hits  the  bottom  plane — a  situa¬ 
tion  never  observed  on  the  Sun. 

One  can  construct  a  rising  flux  rope  with  dipped  apex  by 
making  it  simultaneously  strongly  unstable  with  respect  to  the 
KI  and  TI.  While  such  a  situation  can  be  realized  in  laboratory 
and  corresponding  numerical  simulations  of  filament  eruptions 
dBellan  &  Hansenl  19981: 1 Arnold  et  al.l2008l).  it  appears  very  spe¬ 
cial  and,  hence,  unlikely  to  occur  under  solar  conditions. 

We  conclude  that  the  writhe  of  erupting  filaments  can  be  reli¬ 
ably  estimated  from  the  apex  rotation  in  the  relevant  height  range 
where  the  nonlocal  writhe  dominates,  unless  the  filament  devel¬ 
ops  a  very  particular  shape  with  multiple  maxima  or  minima. 

5.  Summary 

The  availability  of  3D  observations  provided  by  the  STEREO 
spacecraft,  combined  with  recently  developed  analytical  expres¬ 
sions,  facilitates  obtaining  the  writhe  of  helical  structures  in  the 
solar  corona.  In  order  to  explore  the  relevance  of  this  quantity 
for  coronal  phenomena,  we  investigated  how  it  relates  to  the 
projected  shape  of  helical  curves  and  measured  its  evolution  in 
numerical  simulations  of  ideal  MHD  flux  rope  instabilities.  Our 
results  and  their  implications  for  stable  and  eruptive  coronal  ob¬ 
jects  can  be  summarized  as  follows. 


(1)  The  relation  between  writhe  and  projected  S  shape  of  a 
curve  with  both  end  points  in  a  plane  is  not  unique.  It  depends  on 
the  height  of  the  curve  and  on  the  presence  or  absence  of  dips. 
Therefore,  in  principle,  care  must  be  taken  when  associating  a 
sign  of  writhe  to  an  observed  S  shape  on  the  Sun.  However,  we 
demonstrated  that  this  ambiguity  does  not  affect  low-lying  fila¬ 
ments,  as  long  as  their  height  remains  below  about  0.4  their  foot- 
point  distance,  or  soft  X-ray  sigmoids.  This  supports  the  estab¬ 
lished  empirical  rule  which  associates  stable  forward  (reverse) 
S  shaped  structures  low  in  the  corona  with  positive  (negative) 
helicity. 

(2)  Kink-unstable  erupting  flux  ropes  transform  a  far  smaller 
fraction  of  their  twist  helicity  into  writhe  helicity  than  often  as¬ 
sumed.  The  writhe  and  the  number  of  turns  of  a  helical  object 
differ  in  general.  In  particular,  a  simple  leg  crossing  of  a  rising 
filament  in  sky  projection  can  not  be  taken  as  evidence  that  the 
writhe  of  the  filament  is  close  to  unity. 

(3)  Confined  flux  rope  eruptions  tend  to  show  stronger  writhe 
at  low  heights  than  ejective  eruptions  (CMEs),  which  acquire 
writhe  over  a  larger  height  range.  This  argues  against  sugges¬ 
tions  that  the  writhing  facilitates  the  rise  of  unstable  flux  ropes 
through  the  overlying  field,  and  it  implies  that  a  significant  ro¬ 
tation  may  occur  in  a  larger  fraction  of  CMEs  than  suggested 
by  the  apex  rotation  of  associated  erupting  filaments  low  in  the 
corona. 

(4)  Erupting  filaments  which  are  S  shaped  already  before  the 
eruption  and  keep  the  sign  of  their  axis  writhe  (which  is  expected 
if  field  of  one  chirality  dominates  the  source  volume  of  the  erup¬ 
tion),  must  reverse  their  S  shape  in  the  course  of  the  rise.  In  flux 
rope  topology,  the  initial  straightening  up  to  the  reversal  repre¬ 
sents  a  conversion  of  writhe  into  twist,  which  cannot  be  caused 
by  the  helical  kink  instability.  However,  the  writhe  conversion 
can  be  a  mechanism  that  triggers  the  instability  after  the  straight¬ 
ening. 

(5)  The  writhe  of  rising  loops  of  simple  shape  (which  do  not 
have  secondary  maxima  or  minima)  can  be  estimated  from  the 
angle  of  rotation  about  the  direction  of  ascent,  once  the  apex 
height  exceeds  the  footpoint  separation  significantly.  This  pro¬ 
vides  a  convenient  means  to  estimate  the  writhe  of  erupting  fil¬ 
aments  which  rise  towards  the  observer,  even  if  no  3D  observa¬ 
tions  are  available. 

We  emphasize  that  for  flux  tubes  where  the  twist  varies 
significantly  with  radius,  the  interplay  between  twist,  writhe, 
and  helicity  becomes  more  difficult  (iLongcope  &  Malanushenkol 
20081).  Usually  an  average  of  the  twist  over  the  cross  section  of  a 
flux  tube  is  of  interest,  while  the  writhe  is  a  property  of  a  single 
field  line,  in  this  context  the  reference  field  line  of  the  twist. 

This  work  shows  that  writhe  is  a  useful  quantity  in  inter¬ 
preting  S  shaped  coronal  structures  and  in  constraining  eruption 
models.  It  can  straightforwardly  be  computed  for  numerical  data 
and  can  often  be  estimated  from  observations.  It  has  a  range  of 
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further  applications  in  solar  physics,  for  example,  the  study  of 
the  dynamo  (lAsgari-Targhi  &  Berged  2009).  of  flux  ropes  ris¬ 
ing  through  the  convection  zone  (e.g..  IlLinton  et  al.ll  19981).  and 
of  magnetic  loops  connecting  different  active  regions,  including 
transequatorial  loops  (e.g  JChen  et  al.ll2007l). 
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